RANK OF STABLY DISSIPATIVE GRAPHS 



PEDRO DUARTE AND TELMO PEIXE 



Abstract. For the class of stably dissipative Lotka-Volterra systems we prove that 
the rank of its defining matrix, which is the dimension of the associated invariant 
foliation, is completely determined by the system's graph. 



1. Introduction 



In his work "Legons sur la Theorie Mathematique de la Lutte pour la Vie" [T7] 
Volterra began the study of differential equations 



where Xi{t) > represents the density of population i in time t and its intrinsic rate of 
growth or decay. Each coefficient aij represents the effect of population j on population 
i. If ttij > this means that population i benefits from population j. A = (aij) is said 
to be the interaction matrix of the system ([T]). This system of differential equations is 
usually referred to as the Lotka-Volterra equations. 

For two-dimensional systems we can completely analyze the dynamics of ([T]) (see, 
for example, [H]). However we are far from being capable of doing the same for higher 
dimensional Lotka-Volterra systems, in spite of some important results [HI UHl El 121] • 

There is a close relation between studying the dynamical properties of a Lotka- 
Volterra system and the algebraic properties of its interaction matrix, and depending 
on that these systems can be classified in tree main classes: cooperative or competitive; 
conservative; and dissipative. 

For the cooperative and competitive systems , overall results were obtained by Smale 
[H] and Hirsch |3lll], among others, for example, Zeeman [THl [191 [20] , Van Den Driess- 
che et al. [1], Hofbauer et al. [6j, Smith [15] and Karakostas et al. [7|. 

Concerning the conservative systems, the initial investigations are attributed to 
Volterra, who also defined the class of dissipative systems [T7] looking for a gener- 
alization of the predator-prey system. Given A G Mat„xn(I^) and a point q G we 
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can write the Lotka-Volterra system ([T]) as 

f = , (2) 

where X^^glx) = x * A{x — q). The symbol denotes point-wise multiplication of 
vectors in M". We say that system ([2]), the matrix A, or the vector field X^^q, are 
dissipative iff there is a positive diagonal matrix D such that Q{x) = x^ADx < for 
every x G M". Notice this condition is equivalent to x^D~^Ax < 0, because 

x'^D~^Ax = {D-^xfAD{D-^x) = Q{D-^x) < . 

When A is dissipative, ^ admits the Lyapunov function 

n , 

7/ \ \^Xi — qiiogXi , , 

^(^) = 2^ J. ' (3) 

i=l * 

which decreases along orbits of XA,q- In fact the derivative of h along orbits of XA,q is 

n 

^ = E !r^^^ - - ^^-^ = - qVD-'A{x - g) < 0. (4) 

Since h function is proper, ^ determines a complete semi-flow 0^^^ : M" >, defined 
for all t > 0. Let F^.g denote the forward limit set of i.e., the set of all accumulation 
points of (pA^g{x) as t — 7- +oo, sometimes referred to as the system's attractor. By La 

Salle's theorem [Sj we know that for dissipative systems TA,q C [h = 0}. 

The notion of stably dissipative is due to Redheffer et al, who in a series of papers 
in the 80's [I2l [13], [HI IHl HD] studied the asymptotic stability of this class of systems, 
using the term stably admissible systems. Redheffer et al. designated by admissible 
the matrices that Volterra had initially classified as dissipative [TTj. We now give the 
precise definition of stably dissipative system. Given a matrix A G Mat„xn(I^) we say 
that another real matrix A G Mat„xn(I^) is a perturbation of A iff 

= ^ aij = 0. 

We say that a given matrix A, XA,q, or (|2]), is stably dissipative iff any sufficiently small 
perturbation A of A is dissipative, i.e., 

3e > : max \ aij — aij\ < e =^ A is dissipative. 
From the interaction matrix A we can construct a graph Ga having as vertices the n 



species {1, . . . , n}. See definition 3.1 below. An edge is drawn connecting the vertices i 
and j whenever aij 7^ or aji 7^ 0. Redheffer et al. [121 [131 [Til |9l|10] have characterized 
the class of stably dissipative systems and its attractor TA,q in terms of the graph Ga- 
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In particular, they describe a simple reduction algorithm, running on the graph Ga-, 
that 'deduces' every restriction of the form T C {x : = gj}, 1 < i < n, that holds 
for every stably dissipative system with interaction graph Ga- To start this algorithm 
they use the following 

Lemma 1.1. Given a stably dissipative matrix A, there is some positive diagonal ma- 
trix D such that 

n 

WiWj = ^ QiiWi = 0, y i = 1, . . . ,n . 

Since by La Salle's theorem 

Ta,, C J X e M!^ : ^ (x, - qi) {xj - qj) = 

I ij=l * 

it follows that ^ C {x : Xi = qt} for every i = 1, . . . ,n such that ajj < 0. A vertex i 
is colored black, •, to state that g C {x : Xj = qi} holds. Similarly, a cross is drawn 
at a vertex i, ©, to state that F^^^ C {x : X^^^(x) = 0}, which means {xj = const} 
is an invariant foliation under (p'^j^^ : Ta^q All other vertices are coloured white, 
o. Before starting their procedure, as an < for all i, they colour in black every 
vertex i G {1, ...,«} such that an < 0, and in white all other vertices. This should 
be interpreted as a collection of statements about the attractor F^i g. The reduction 
procedure consists of the following rules, corresponding to valid inference rules: 

(a) If j is a • or ©-vertex and all of its neighbours are •, except for one vertex /, 
then colour / as •; 

(b) If j is a • or ©-vertex and all of its neighbours are • or ©, except for one vertex 

then draw © at the vertex /; 

(c) If j is a o-vertex and all of is neighbours are • or ©, then draw © at the vertex 
j- 

Redheffer et al. define the reduced graph of the system, 01{A), as the graph obtained 
from Ga by successive applications of the reduction rules (a), (b) and (c) until they 
can no longer be applied. In [11] Redheffer and Walter proved the following result, 
which in a sense states that the previous algorithm on Ga can not be improved. 

Theorem 1.2. Given a stably dissipative matrix A, 

(a) // "JllA) has only •-vertices then A is nonsingular, the stationary point q is 
unique and every solution of converges, as t oo, to q. 

(b) // !3i(y4) has only • and (B-vertices, but not all •, then A is singular, the sta- 
tionary point q is not unique, and every solution of has a limit, as t oo, 
that depends on the initial condition. 
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(c) //3^(y4) has at least one o-vertex then there exists a stably dissipative matrix A, 
with = Ga, such that the system ^) associated with A has a nonconstant 
periodic solution. 

In a very recent paper, Zhao and Luo [21] gave necessary and sufficient conditions 



for a matrix to be stably dissipative, see proposition 3.5 



The previous theorem tells us that when Jl{A) has only w-vertices, then the matrix 
A has always maximal rank, rank(74) = n. In case Jl{A) has only • and ©-vertices, 
then rank of the matrix A equals the dimension of an invariant foliation. This led us 
to establish that the rank of any stably dissipative matrix only depends on its graph. 
In particular the same is true for the dimension of the invariant foliation of any stably 
dissipative system. See theorem |4.2| and corollary |4.3| of section |4j 



2. Dissipative Systems and Invariant Foliations 

Assume the Lotka-Volterra field XA,q{x) = x* A{x — q), defined in ([2]), is associated 
with a dissipative matrix A with rank{A) = k. Let W he a [n — k) x n matrix whose 
rows form a basis of 

Ker(y4^) = { a; G M" : x'^A = } . 

Define the map g -.Rl ^ M"~^ g{x) = Wlogx, where logx = (logxi, . . . , logx„). 
Denoting by = diag(xi, . . . , x„) the diagonal matrix, the jacobian matrix of g, 
Dgx = W D~^, has maximal rank n — k. Hence g : M" — )■ M""'^ is a submersion. 

We denote by 3^ a the pre-image foliation, whose leaves are the the pre-images g~^{c) 
of g. By a classical theorem on Differential Geometry, each non-empty pre-image g~^{c) 
is a submanifold of dimension k. Recall that the dimension of a foliation is the common 
dimension of its leaves. A foliation 5" is said to be invariant under a vector field X if 
X{x) G TxJ-' for every x, where T^J^ denotes the tangent space at x to the unique leaf 
of J-" through X. This condition is equivalent to say that the flow of X preserves the 
leaves of J-". 

Proposition 2.1. // A is dissipative then the foliation Ta is XA,q-invariant with 
dim(Jvi) = rank{A). 
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Proof. We have 

Dg^{XA,q{x)) = Dg^{D^A{x - q)) 
= WD;^D,A{x - q) 
= WA(x - = , 

because 1^74 = 0. Hence XA,q{x) e T^Ta and Ta is X^,q-invariant. □ 
Defining the symmetric and skew- symmetric parts of a matrix A by 

^sym ^ + ^ ^skew ^ 

2 2 
the following decompositions hold 

A = A'^"^ + A'^"""" and A^ = A'^"^ - A'^'''" . 

The following lemma is a simple but key observation 
Lemma 2.2. Given a dissipative matrix A, 

Ker(^) = Ker(^^) = Ker(A^y'^) n Ker(A^^"^). 



Proof. It is obvious that Ker(^) D Ker(^"2^"^) n Ker(^"^'="'). On the other hand, if 

V e Ker(A) then Av = 0, and hence v'^A^^^'v = v'^Av = 0. Since A^^"" < 0, it follows 
that V e Ker(A^2/m) Observing that A'''^'" ^ A- A'^^ , we have v e Ker(A^*^^^). Thus 

V e KeiiA^y"") n Ker(A"^^'^). We have proved that Ker(A) = Ker(^^2/m) p Ker(A^'=^^). 
Analogously, Ker(^^) = Ker(^*2/m-) ^ Ker(^"'=^"'). □ 



Define 

EA,q = {xeWl : A{x-q) = Q} (5) 
to be the affine space of equilibrium points of Xa^q- 



Theorem 2.3. Given a dissipative matrix A, each leaf of J^a intersects transversely 
EA,q in a single point. 
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Proof. Let V he a kxn matrix whose rows form a basis of Ker(y4)-'-, the space generated 
by the rows of A. With this notion, A{x — g) = is equivalent to V{x — q) = 0. Let 
W be the {n — k) x n matrix used in the definition of J-'a, whose rows form a basis 

of Ker(y4"^). With W and V we form the n x n matrix U = ^ 



V 



. By lemma 



2.2 



Ker(A-^) = Kei{A), and this implies that U is nonsingular. The intersection of a lea: 
g'^{c) of J-'a with the equlibria set EA,q is described by the system 



X e g (c) n EA,q ^ 



W log X = c 
V{x -q) = 



Substituting u = logx, this system becomes equivalent to 

Wu = c 



y{e - g) = ■ 

In order to see that each leaf of Ta intersects the equilibria set EA^q at a single point, 
it is enough to see that 

Wu = c { Wu' = c 

Vie'' -q) = j V(e«' - q) = 

imply u = u'. By the mean value theorem for every i G {1, 
Ui e such that 

which in vector notation is to say that 

e"" - e"' = De4u - u') = e" * (m - u'). 



n} there is some 



Hence 



W{u -u')=0 
Vie - e"') = 

w 



W{u - m') = 
VDAu-u!\ = 



VD^u 



[u — u 



^ U 



' I 











{u-u') = . 



Therefore, as 














is non-singular, we must have u = u' . The transversality 



follows from the nonsingularity of U 



' I 











□ 
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3. Stably Dissipative Systems 

Definition 3.1. Given a dissipative matrix A we define its graph as Ga = (^Aj-A-a), 
with Va = V,UVo, V, = {1 <i <n : an < 0} and Vo = {I < i < n : a^j = }. A 

pair {i,j), with i ^ j , is an edge in Aa iff aij or aji ^ 0. 



Consider a simple graph G = (V, ^1) whose vertices are colored in black and white, 
meaning there is a decomposition V = V, U Vo of the vertex set. We shall say that 
a vertex in Vo is a o-vertex, while a vertex in V, is a •-vertex. Such graphs will be 
referred as black and white graphs. 

In [H] Redheffer and Walter gave the following important property of stably dissi- 
pative matrices in terms of their associated graph. 



Lemma 3.2. If A is a stably dissipative matrix then every cycle of Ga has at least one 
strong link (•—•)• 

Proof. The proof is by contradiction, for otherwise we could perturb A into a non 
dissipative matrix. □ 



Definition 3.3. A black and white graph G is called stably dissipative iff every cycle 
of G contains at least a strong link, i.e., an edge between •-vertices (•—•)• 

The name 'stably dissipative' stems from the use we shall make of this class of graphs 
to characterize stably dissipative matrices. See Proposition |3. 5 below. 

Let us say that a dissipative matrix A G Mat„xn(I^) is almost skew- symmetric 
iff ttij = —aji whenever an = or ajj = 0, and the quadratic form Q{xk)kev. = 
jev. '^ij ^3 is negative definite. 



Proposition 3.4. Given a dissipative matrix A G Matnxn(IR); there is a positive di- 
agonal matrix D such that a^j dj = —aji di whenever an = or ajj = 0, and for every 
Xfc G M with k eV,, Xlijev. CLijdiXiXj < 0. 

Proof. Let be a positive diagonal matrix such that for all x G M", Q{x) = ^ . j aijdiXiXj < 
0. Assuming an = 0, choose a vector x G M" with Xi = and Xk = for every k ^ i,j. 
Then for every Xj G M, 

which implies that atjdj + ajidi = 0, and everything else follows. □ 
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Recently, Zhao and Luo |2T] gave a complete characterization of stably dissipative 
matrices. 



Proposition 3.5. Given A G Mat„xn(I^); A is stably dissipative iff Ga is a stably 
dissipative graph and there is a positive diagonal matrix D such that AD is almost 
skew- symmetric. 



Proof. We outline the proof. Assuming A is stably dissipative, by lemma 



^ Ga 
which 
AD 



3.4 



is stably dissipative. Take a diagonal matrix D > according to lemma L_ 
implies that Q{xk)k£V. = J2ijev. '^ijdiXiXj is negative definite. By proposition 
is almost skew-symmetric. 

Conversely, assume Ga is stably dissipative, assume AD is almost skew-symmetric, 
and take A = {(iij) some close enough perturbation of A. Let Ga be the partial graph 
of Ga obtained by removing every strong link (•—•). Because A is stably dissipative, 
the graph Ga has no cycles. We partition the vertex set {1, . . . ,n} as follows: Let Vq 
be a set with an endpoint in each connected component of Ga- Recursively, define Vk 
to be the set of all vertices such that there is an edge of Ga connecting it to a vertex 
in Vk-i- By construction, we have a map i i— ?■ i\ associated with this partition, such 
that i' G Vfc_i for every z G with > 1, and Ga has an edge connecting i with i' . 
Then we consider the diagonal matrix D = diag{dj) whose coefficients are recursively 
defined hy di = di if i E Vq, and 



di = —di' for i G Vfc with A; > 1 . 



di'i 

It follows by induction in k that di > for every i G Vfc. For k = this is automatic. 
Assuming this holds in Vk-i, for any i G we have d^r > 0, which implies that di > 
because da' and di'i have opposite signs. The diagonal matrix D is close to D because 
A is near A. Therefore, by continuity, the quadratic form Q{xk)kev. = J2ijev. dijdiXiXj 
is negative definite. On the other hand, by definition of di, di dai + dii di'i = 0. Since 
every o-o and o-« connection links i with i' for some i ^ Vq, this implies that for every 
(xj) G M", J2'ij=i (^ijdiXiXj < 0. Hence A is dissipative, which proves that A is stably 
dissipative. □ 
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4. Main Results 

Definition 4.1. Given a stably dissipative graph G, we denote by SD(G) the set of all 

dissipative matrices A with Ga = G. 

Our main theorem is 

Theorem 4.2. Given a stably dissipative graph G, every matrix A G SD(G') has the 
same rank. 

By this theorem we can define the rank of a stably dissipative graph G, denoted here- 



after by rank(G), as the rank of any matrix in SD(G). Together with proposition 2.1 
this imphes that 

Corollary 4.3. Given a stably dissipative graph G, for every matrix A G SD(G), any 
stably dissipative Lotka-Volterra system with matrix A has an invariant foliation of 
dimension = rank(G). 

Definition 4.4. We shall say that a graph G has constant rank iff every matrix 
A G SD(G) has the same rank. 



With this terminology, theorem 4^ just states that every stably dissipative graph 
has constant rank. 



5. Simplified Reduction Algorithm 

As before, let TA,q denote the forward limit set of (|2]), i.e., the set of all accumu- 
lation points of (pA,qi^) as t — > +oo. We say that a species i G {1, . . . ,n} is strongly 
dissipative iff TA,q C {x G M" : Xi = qi }, or equivalently limf_>.+oo = 'iii fo'^ 

all X G M". Similarly, we say that a species i G {1, . . . , n} is weakly dissipative iff 
limt_5.+oo 0A,g,j(a;) exists, for all x G M" . With this terminology, the algorithm of Red- 
heffer et ai, described in the introduction, is about the determination of all 'strongly' 
and 'weakly dissipative' species of the stably dissipative system Since the algo- 
rithm runs on the graph Ga, the conclusions drawn from the reduction procedure hold 
for all stably dissipative systems that share the same graph Ga- 



The following proposition is a slight improvement on item (b) of theorem 1.2 



Proposition 5.1. IfJllA) has only • and (B-vertices then the system has an invariant 
foliation with a single globally attractive stationary point in each leaf. 
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Proof. Combine theorem |1.2[ b) with proposition 2.3 □ 



In fn\ Redheffer and Zhiming make the following remark: 

Remark 5.2. Let A be dissipative and let every vertex o in Ga be replaced arbitrarily 
by ©. Then A is nonsingular iff, by algebraic manipulations, every vertex can then 
be replaced by •. 

We shall explain this remark in terms of a simpler reduction algorithm. Denote 
by EA,q the set of all equilibria of ([2]), EA,q = {x E M" : A(x — q) =0}. Let us 
say that a species i G {1, . . . ,n} is a restriction to the equilibria of XA,q whenever 
EA,q C {x G M" : Xi = qi}. Notice that every strongly dissipative species is also a 
restriction to the equilibria of ^a,?- Think of colouring i as black as the statement 
that i is a restriction to the equilibria of XA,q- Notice that at the begining of the 
reduction algorithm, described in the introduction, the weaker interpretation that all 
black vertices correspond to restrictions to the equilibria is also valid. If we simply 
do not write ©-vertices, but consider every o-vertex as a ©-vertex, then the reduction 
rules (b) and (c) can be discarded, while the first rule, (a), becomes 

(R) If all neighbours of a vertex j are ^-vertices, except for one vertex k, then we 
can color k •-vertex. 

The idea implicit in the remark above is that (R) is a valid inference rule for the 
weaker interpretation of the colouring statements above. Assuming that every o-vertex 
is a ©-vertex amounts to looking for restrictions on the equilibria set EA,q instead of the 
attractor g. Let us still call reduced graph to the graph, denoted by ^R^(G), obtained 
from G by successively applying rule (R) alone until it can no longer be applied. The 
previous considerations show that 

Proposition 5.3. Given a stably dissipative matrix A, every •-vertex of'Jl^:{GA) is a 
restriction to the equilibria of XA^q- 

We shall write "Ji^iG) = {•} to express that all vertices of 'Ji^:{G) are •- vertices. 



Corollary 5.4. If G is a stably dissipative graph such that = 
matrix A G SD(G') is nonsingular. In particular G has constant rank. 



{•} then every 



Proof. Given A G SD(G'), by Proposition 5.3 we have EA,q 
implies that A is nonsingular. 



{q}, which automatically 

□ 



In fact, the converse statement of this corollary holds by remark 



5.2 
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Proposition 5.5. Let A be stably dissipative matrix. If A is nonsingular then '3i^,{G a] 
{•}■ 



6. Proofs 

We call any extreme o- vertex of G a o-endpoint of G. 

Lemma 6.1. // a stably dissipative graph G has no o-endpoints then I3^*(G) = {•}. 

Proof. Let G be a stably dissipative graph with no o-endpoints. Assume, by contradic- 
tion, that 7^ {•}. We shall construct a cycle in Jl^:{G) with no •— • edges. Since 
every o- vertex of is also a o- vertex of G, this will contradict the assumption 
that G is stably dissipative. In the following construction we always refer to the vertex 
coloring of Jl^^lG). Take jo to be any o- vertex. Then, given jk take a neighbouring 
vertex jk+i to be another o-vertex, if possible, or a w-vertex otherwise. While the path 
is simple (no vertex repetitions) it can not end at some o-endpoint, and it can not 
contain any •— • edge because whenever we arrive to a w-vertex from a o-one we can al- 
ways escape to another o-vertex. In fact, no •- vertex can be linked to a single o-vertex 
since otherwise we could reduce it to a •- vertex by applying rule (R). By finiteness this 
recursively defined path must eventually close, hence producing a cycle with no •— • 
edges. □ 



Given a stably dissipative graph G and some o-endpoint z G Vo, we define the 
trimmed graph Ti{G) as follows: Let G V be the unique vertex connected to i by 
some edge of G. Then Ti{G) is the partial graph obtained from G by removing every 
edge incident with i'. The trimming operation preserves the stable dissipativeness of 



graph, i.e., Ti(G) is stably dissipative whenever G is. This follows by definition 3.3 
because Ti{G) is obtained by removing some edges from G. 




G Ti[G) 



Figure 1. A graph G and it's trimmed graph Ti{G). 
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Similarly we define the trimmed matrix Ti{A) as follows: Annihilate every entry of 
row i', except for Oj/j and 0,/,/, and annihilate every entry of column i', except for an/ 
and ai'i'. 



A 



■ ■ * • 
an' ■ 

■ ■ * ■ 

* ai'i * o^i'i' * 

■ ■ * ■ 



■ • • 

an' ■ 

■ ■ ■ 

ai'i ai'i' 

■ ■ ■ 



The above represent entries of A that are annihilated in Ti[A). 



Lemma 6.2. Let i G Vq be some o-endpoint of a stably dissipative graph G. 
IfAeSD{G) then Ti{A) e SD{Ti{G)) and rank(Ti(A)) = rank(A). 

Proof. Take A G SD(G) and let A' = Ti{A), where i is some o-endpoint. Denote, 
respectively, by colj and rowj the j^^ column and the j^^ row of A, and denote by col^- 
and Tow'j the j*^ column and the j^^ row of the trimmed matrix. Since i is a o-endpoint, 
an' is the only nonzero entry in rowj, and Oj^j is the only nonzero entry in colj. Then 
the trimmed matrix A' is obtained from A by applying the following Gauss elimination 
rules, either simultaneously or in some arbitrary order 



rowl- := row, ^row,- 



dii' 



col'- := coL- — coL j ^ i' . 

' ' ai',i 

Because Gauss elimination preserves the matrix rank we have rank(y4') = rank(y4). 
To finish the proof, it is enough to see now that A' is stably dissipative. We use 



proposition 3.5 for this purpose. First, Ga' = Ti{G) is stably dissipative as observed 



above. Let D be a positive diagonal matrix such that AD is almost skew-symmetric. In 



view of proposition 3.5, we only need to prove that A' D is also almost skew-symmetric. 



Notice that G and Ti{G) share the same black and white vertices. If a'^^ = or 
a'jj = then also akk = or ajj = 0. Hence, because AD is almost skew-symmetric, 
akj dj = —ajk dk. Looking at the Gauss elimination rules above, we have a'j^j = akj and 
^'jk — '^jky or else a'f^j = a'-j^ = 0. In either case we have a'f^^dj = —a'j^^dk. Finally, we 
need to see that Q'{xe)£ev. = J2kjev. ^'kj^j^kXj is a negative definite quadratic form. 
If i' is a o- vertex then Q'{xi)ifzv, = J2kj&v. ^kjdj XkXj is negative definite because AD 
is almost skew-symmetric. Otherwise, if i' is a w-vertex, given a nonzero vector {xi)£^\i. 
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we define (x^)^^^. letting x'^ = xe for £ ^ i' , wliile x[, = 0. Tlien 

Q'{xi)i(zv. = tti'ii di'xl + ^ akjdjx'^x'j < , 

^ V ' 

=Q(4)f6V.<0 

since {xe)eev. 7^ implies that either Xj/ 7^ or else (x^)^^^. 7^ 0. This proves that Q' 
is negative definite. □ 



As a simple corollary of the previous lemma we obtain 

Lemma 6.3 (Trimming lemma). Let i G Vo be some o-endpoint of a stably dissipative 
graph G. If Ti{G) has constant rank then so has G, and rank(G) = rank(Tj(G)). 



We can now prove theorem 4.2 



Proof Define recursively a sequence of graphs Gq, Gi, . . . , Gm, with Go = G, and where 
Gj+i = Tj.{Gi) for some o-endpoint ji of Gj. This sequence will end at some graph Gm 
with no o-endpoint. By Lemma 6.1 we have 3^^,(Gm) = {•}• The connected components 
of Gm are either reducible to •- vertices by iteration of rule (R), or else composed by 
o-vertices alone. Since the o-components can not be trimmed anymore, they must be 
either formed of a single o- vertex, or else a single 0—0 edge. By corollary |5.4[ Gm has 



constant rank. Finally, applying inductively Lemma 6^ we see that all graphs Gj have 
constant rank. Hence G, in particular, has constant rank. □ 



Remark 6.4. The previous proof gives a simple recipe to compute the rank of a graph. 
Trim G while possible. In the end, discard the single o-vertex components and count 
the remaining vertices. 
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Table 1 . Some graph trimming examples. 
Original graph Trimmed graph Graph rank 



o o 



7. Trimming Effect on Dynamics 

In this last section we use an example to describe the effect of trimming a stably 
dissipative matrix on the underlying dynamics. 
Consider the system 



E 



Xi 


= Xi{{x2 - 1) + 


(X7 - 1)) 


X2 


= X2{-2{xi - 1 


) + (^3 - 1)) 


X3 


= -X3{X2 - 1) 




X4 


= X4{{X5 - 1) - 


(^7 - 1)) 


X5 


= a;5(-2(,T4 - 1 


) + (^6 - 1)) 


Xq 


= -X6{X5 - 1) 




X7 


= X7{-{xi - 1) 


+ {X4 - 1) - 



with interaction matrix 



A 






1 














1 


-2 





1 

















-1 





























1 





-1 











-2 





1 

















-1 








-1 








1 








-1 
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and fixed point g=(l,l,l,l,l,l,l,). E has associated graph Ga represented in figure 

m 



o — o — o — • — o — o — o 

3 2 1 7 4 5 6 



Figure 2. Associated graph of matrix A, G{A). 




Figure 3. Phase portrait of a system E. 

The null space, Ker(y4), is generated by the vector (1,0,2,1,0,2,0). Hence the 
foliation J-", with leaves J-'c given by 

J'c = {x EM7 : log Xi + 2 log X3 + log X4 + 2 log Xq = c} , 

is an invariant foliation with dimension 6 in M7 . The system's phase portrait is repre- 
sented in figure |3| being the atractor a 3-plan transversal to J-" given by 

r = {s G M'' : Xi = X4, X2 = X5, Xs = Xq, X7 = 1} . 

The intersection of each leaf J-'c with F is a surface Sc given by 

5c = -Fc n F = {{xi,X2,X3,xi,X2,X3, 1) : loga;i + 21ogX3 = -}, 

wich is foliated into invariant curves by the level sets of h, defined in ([3]). Note that Sc 
corresponds to an invariant leaf of the conservative system with graph o — o — o. 

With the first trim on G we get the graph Tq{G) represented in figure |4j 
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o — o — o — • — o o — o 

3 2 1 7 4 S 8 

Figure 4. The trimmed graph of G, Tq{G). 

This corresponds to annihilate the entries (4,5) and (5,4) of the original matrix 
A. Notice that the components and Xq of the system are independent of the rest. 
Hence the dynamics of this new system is the product of two independent LV systems 
represented in figure [5] 




Figure 5. Representation of the dynamic of the system Ei. 

The five dimensional system on the left of figure [5] has a straight line of equilib- 
ria. Moreover it leaves invariant a foliation of dimension four with a single globally 
attractive fixed point on each leaf. The right-hand side system is a typical conservative 
predator-prey. 

Now we have two different possibilities of trimming the graph Tg(G'): we can ei- 
ther choose the o-endpoint 3 or else 4. In the first case we get the graph T^{Tq{G)) 
represented in figure [6} whose dynamics is illustrated in figure [7} 

o — o o — • — o o — o 

3 2 1 7 4 5 6 

Figure 6. The trimmed graph r3(r6(G)) of Tq{G). 

The three dimensional system in the middle of figure [7] has a straight line of equi- 
libria. Moreover it leaves invariant a foliation of dimension two with a single globally 
attractive fixed point on each leaf. The left and right-hand side systems are typical 
conservative predator-preys. 

In the second case we get the graph T^iTQiG)) represented in figure |8| whose dy- 
namics is depicted in figure [9] 
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Figure 7. Representation of the system's dynamics associated to the graph T^(Tq{G)) . 



o — o — o • — o o — o 

3 2 1 T 4 5 8 



Figure 8. The trimmed graph r4(r6(G)) of r6(G). 
















X 






X 

















Figure 9. Representation of the system's dynamics associated to the graph T4^{Tq{G)) . 

Here, the left-hand side three dimensional system is conservative, leaving invariant 
a foliation of dimension two transversal to a straight line of equilibria. The middle and 
right-hand side systems are typical predator-prey, respectively, dissipative and conser- 
vative. 



Trimming T4{Tq{G)) choosing the o-endpoint 1 we get the graph Ti{T4(Tq{G))) rep- 
resented in figure [10} whose dynamics is a product of three predator-prey systems, 
illustrated in figure [TTl with a one dimensional system consisting of equilibria. 



Notice that by trimming T3(Tq{G)) we obtain an isomorphic graph to the one in 



figure 10 



o o — o 

3 2 1 



o o — o 

4 5 8 



Figure 10. The trimmed graph Ti{Ti{Te{G))) of Ti{TeiG)). 
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Figure 11. Representation of the system's dynamics associated to the graph Ti{T4,{Tq{G))) . 
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